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We propose mutant versions of quantum mechanics constructed on vector spaces over the finite 
Galois fields GF(3) and GF(9). The mutation we consider here is distinct from what we proposed 
in previous papers on Galois field quantum mechanics. In this new mutation, the canonical expres- 
sion for expectation values is retained instead of that for probabilities. In fact, probabilities are 
indeterminate. Furthermore, it is shown that the mutant quantum mechanics over the finite field 
GF(9) exhibits super-quantum correlations (i.e. the Bell-Clauser-Horne-Shimony-Holt bound is 4). 
We comment on the fundamental physical importance of these results in the context of quantum 
gravity. 
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I. INTRODUCTION 

Quantum descriptions of physical systems begin with 
the introduction of a vector space, generally defined over 
the complex number field C, with elements of the space 
associated with states of the physical system under con- 
sideration. This space, in the traditional approach, is 
assumed to be a Hilbert space "H, which for N- level sys- 
tems is H — C N . The Hilbert space rl possesses a natural 
inner product WxH^C, which we denote 



(|a),|/3» e C, \a),\(3)eH. 



(1) 



It is customary to associate a dual- vector (a\ G H* to 
each vector \a) <G rl, with the same label a, via 



so that 



(a\ = (\a), ) , 



(a\(3) = (\a),\/3)) 



(2) 



(3) 



The presence of the inner product allows the definition 
of hermitian conjugation of linear operators via 



(|a),i|/3» = (i^a),!/?)) , 



(4) 



and that of hermitian operators, to which physical ob- 
servables are associated, via A^ = A. It also allows for 
the definition of unitary operators via 



(U\a),U\l3)) = (|a>,|/?» 



(5) 



under which the states are assumed to evolve. 

There are two ways for the quantum description to 
make contact with physical reality. In the first approach, 
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possible outcomes of a measurement of an observable A 
are assumed to be given by its eigenvalues. Let us denote 
the eigenvector associated with eigenvalue a by \a): 



(A-a)\a) = . 



(6) 



When the system is in the state represented by \ip) E H, 
the probability of obtaining the outcome a € M as a 
result of a measurement of A is given by 



P(a\ 



(7) 



where the sum in the denominator runs over all the eigen- 
values of A. The hermiticity of A ensures that its eigen- 
values are all real, and that the eigenvectors are mutually 
orthogonal and complete. Normalizing the state vector 
and the eigenvectors of A so that 



1 



5 a p , 



(8) 



the above expression for the probability reduces to 

p(c#)HHV>| 2 . 

There is an alternative way of making contact with 
reality, which is equivalent to the one above for conven- 
tional treatments. One begins with the quantity 



WW 



(9) 



which is real for hermitan A, and interprets the result 
as the expectation value for the associated observable 
in the state \tp). If (ip\ip) = 1, the expression reduces 
to (V'l^lV')- Note that this is the standard approach in 
Quantum Field Theory (QFT) where all physical predic- 
tions are expressed in terms of A^-point correlation func- 
tions, i.e., the vacuum expectation values of the products 
of N field operators. This is most explicit in the path in- 
tegral formulation of QFT. To recover the probabilistic 
interpretation of the first approach, one asserts that the 
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probability for obtaining the outcome a for the measure- 
ment of A on the state \4>) is given by 



P(a\ 



\S(A-a)\ 



(10) 



No absolute values are invoked, and attention is shifted 
to moments of the relevant observable operator in the 
state in question; in particular we do need the expec- 
tation values of powers of the operator. For canonical 
quantum descriptions using the Hilbert space H, these 
two starting points lead to identical results. 

The situation however changes when the underlying 
space is not a Hilbert space. Indeed, for spaces for which 
the inner product is ill-defined, one can expect different 
outcomes for these two approaches. 

In Refs. [1] and [2] (inspired by [3], [4], and [5]), we 
have explored the possibility of discretizing the fields over 
which the vector space is defined but retaining the physi- 
cal interpretation provided by the first approach, namely, 
the definition of probabilities via Eq. (7). The fields we 
considered were finite Galois fields GF(p n ), where neN 
and p is a prime number. For the n = 1 case, they are 
GF(p) = Z p = Z/pZ. Vector spaces over GF(p n ) do 
not have inner products since GF(p n ) is not an ordered 
field 1 , preventing any bilinear map to GF{p n ) from being 
positive-definite (or non-negative) in a natural way. 

However, it was recognized that for Eq. (7) to make 
sense, the dual- vectors that appear in the expression only 
need to constitute a basis for the dual- vector space with a 
possible outcome of a measurement associated with each 
one. The usual pairing of dual- vectors with vectors via 
the inner product is inessential. Indeed, all the inner 
product does, in a sense, is connect the two approaches 
via the property 



= (|V), |a» = (\a), 



(ii) 



so that we can write, 



Eq"|W)| = £qW)*<*W) 

£ Q (V>|a)aHV> (V|%> 



E^IWIV) " (V#> ' 



where we have made the identification 
A = ^ |a)a(o:| . 



(12) 



(13) 



Thus, for the first approach, inner products are not nec- 
essary, and once a basis of the dual- vector space and the 
associated set of outcomes is specified, we have an 'ob- 
servable.' 



1 Ordered fields are fields on which on ordering can be imposed 
that respects both addition and multiplication. 



To make contact with the outcome of measurements 
and probability distributions, we need a map from the 
Galois field to that of non-negative reals. It is essential 
that this map preserves products, which is necessary to 
distinguish entangled states from product ones, and also 
for the actions of symmetry groups on the Galois field. 
This is achieved in [1] and [2] through an absolute value 
function. Eq. (7) can be used as is to define the prob- 
ability of each outcome via the absolute value function 
from GF{p n ) to R given by 



if k = , 
if k ^ . 



(14) 



Here, numbers and symbols with underlines are used to 
denote elements of GF(p n ), to distinguish them from el- 
ements of R. Note that this function is product preserv- 
ing, i.e. \k£\ — \k\\£\, which is essential for probabilities 
of product states to factorize. Applying this formalism 
to 2-level systems, we constructed spin-like observables 
for which the measurement outcomes were ±1 £ R, and 
calculated the Clauser-Horne-Shimony-Holt (CHSH) [6] 
(see also [7-10]) bound for the model and found that it 
was two, despite the fact that no hidden variable mimic 
could reproduce the model's predictions. For details, see 
Refs. [1] and [2]. 

In this paper, we explore consequences of starting with 
the second approach to interpretation, namely, the defi- 
nition of expectation values via Eq. (9). Again, we con- 
sider vector spaces over the finite Galois field GF(p n ), 
which do not have inner products. Thus, the concepts 
of normalizability of states, hermiticity of operators, and 
a dual-vector as a hermitian conjugate of a vector, must 
all be reexamined before we can apply Eq. (9). Further- 
more, working in a vector space over GF(p n ), the expres- 
sion (-01 A\ip) = (row vector) • (matrix) • (column vector) 
will generically lead to an element of GF(p n ), which must 
be mapped to an element of R if the result is to represent 
the expectation value of a measurement of a physical ob- 
servable. While we obtain results similar to our earlier 
ones for certain fields, we discover significant differences 
in others. 

In the following, we will address these points one by 
one and define a 'mutant' QM on vector spaces over the 
fields GF(3) = Z 3 and then GF(9) = Z 3 [i], where i 
is the solution to the equation x 2 + 1 = 0, which is 
irreducible in GF(3). In both cases, we will find that 
(ip\A\il)) £ GF(3) by construction, which will be mapped 
to a number in R. Because we are looking at the expecta- 
tion values of observables, the range of this map need not 
be restricted to the non-negative reals as in the case of 
the absolute value function. In Appendix B we show that 
the requirement that this map preserve products and ac- 
tions of symmetry groups determines the map uniquely. 
It is the use of this map for specific expectation values, 
instead of the absolute value function on brackets, that 
distinguishes between the two approaches to interpreta- 
tion. We will show below that the connection to proba- 
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bilitics given by Eq. (10) for canonical QM is no longer 
valid. In fact, individual probability distributions are not 
fixed in our approach, giving rise to indeterminacies be- 
yond those of canonical QM. Our earlier result in Refs. [1] 
and [2] that the CHSH bound for spin-like systems over 
Galois fields cannot be larger than 2 was predicated upon 
using the first approach starting with Eq. (7). We will 
find that in the second approach, the CHSH bound for 
the GF(3) case is also 2. For the GF(9) case, however, 
the CHSH bound is 4, the maximum possible value. As 
far as we are aware of, this is one of the first explicit 
examples of a non-trivial super-quantum theory. 

Before we proceed to the heart of the matter, we note 
that consideration of discrete mathematical structures is 
not only relevant from an academic point of view. We 
note that such considerations have been seriously under- 
taken in various approaches to the quantum structure of 
space and time, i.e. in various forms of quantum gravity. 
The more complete literature can be found in [11]. 

The outline of the paper is as follows: In section II we 
introduce what we call biorthogonal quantum mechan- 
ics, and in section III we present a few examples of this 
construction. Then in section IV we consider the CHCH 
bound and find an explicit example of a super-quantum 
theory. In section V we show that in such a theory prob- 
abilities are indeterminate. We close in section VI with 
detailed comments about the physical relevance of our 
results. Various details not covered in the main text arc 
presented in two appendices. 



II. BIORTHOGONAL QUANTUM MECHANICS 

In order to adopt the definition of expectation val- 
ues via Eq. (9) onto a vector space over the Galois field 
GF(p n ), one must define the analogue of hermitian con- 
jugation of vectors and linear operators without reference 
to an inner product. In this section, we demonstrate that 
this can be accomplished via biorthogonal systems [12]. 

In the following, we restrict our attention to the Galois 
fields GF(p n ) with p = 3 mod 4 and n = 1 or 2. As we 
will see below, this restriction allows our formalism to 
maintain a close parallel to quantum mechanics defined 
on vector spaces over R (n = 1 case) or C (n = 2 case). 

A. Biorthogonal Systems 

As in the previous section, elements of the finite Ga- 
lois field GF{p n ) are denoted by underlined symbols and 
numbers to distinguish them from elements of R or C. 
The A^-dimensional vector space over GF(p n ) is denoted 
V(N,p n ). A biorthogonal system is a set consisting of 
a basis j|l), |2), • • • , \N)} of the vector space V(N,p n ), 
and a basis {(1|, (2|, • • • , (N\} of the dual vector space 
V(N,p n )* such that 

(r\s)=6 rs , r,* = l,2,--. ,N, (15) 



where 

5™ = (? * r *'' (16) 

II it r = s . 
Such a system can be constructed as follows. 

B. Dot Product 

First, denoting the fe-th element of the vector \a) € 
V(N,p n ) as a k £ GF(p n ), define the 'dot product' in 
V(N,p n ) as 

N 

\a)-\b) = Y.€h e GF(p n ) . (17) 
k=i 

Raising an element to the p-th power is semilinear in 
GF(p n ) since 

(a + bf = {a p + b p ) (18) 

in a field of characteristic p. When n — 1, it is an identity 
transformation due to Fermat's little theorem 

aP- 1 = 1 mod p , Va G Z . (19) 

For the case n — 2, p = 3 mod 4, it is an analogue of 
complex conjugation in C. To see this, first note that the 
equation 

x 2 + 1 = (20) 

is irreducible in GF(p) = 7L p if p = 3 mod 4. 2 Denote the 
solutions to this equation as ±i. Adjoining i to GF(p) — 
1 p gives us GF(p 2 ) = Z p [i]. Elements of this field can 
be expressed as a + i b, where a, b € Z p . Then 

(a + ibf = a p + i p b p = a-ib. (21) 

Furthermore, 

(a + ib) p (c + id) = (ac + bd) + i(ad - be) , 

(c + idy (a + ib) = {ac + bd) - i(ad - be) , (22) 

in particular, 

(a + ib) p (a + ib) = a 2 + b 2 £ Z p . (23) 

Therefore, an d |&H a ) are 'complex conjugates' of 

each other, while \a}\a) is 'real.' Thus, whenp = 3mod4, 
the fields GF{p) = 7L p and GF(p 2 ) = Z p [i] take on the 
roles of K and C. 

In the following, when we say GF(p n ), we will mean 
either GF(p) or GF(p 2 ) with p — 3 mod 4 unless stated 
otherwise. Also, borrowing from standard terminology, 
we will say that two vectors in V(N,p n ) are 'orthogonal' 
to each other when they have a zero dot product, and 
that a vector is 'self-orthogonal' when it is orthogonal to 
itself. 



2 x 2 + 1 = is reducible for p = 2 or p = 1 mod 4 since in those 
cases p — 1 will be a solution. 



4 



C. Conjugation of Vectors 



E. Physical States 



Next, choose a basis {|1),|2),- 
such that: 



\N)} for V(N,p n ) 



\r)-\s) 




if r = s , 
if r 7^ s , 



(24) 



that is, all the basis vectors are orthogonal to each other, 
but none are self-orthogonal. Let us call such a basis 
an 'ortho-nondegenerate' basis. The simplest example 
of an ortho-nondegenerate basis would be such that the 
r-th element of the s-th vector is given by S rs , proving 
that such a basis always exists. On the other hand, not 
all bases satisfy this condition since V(N,p n ) typically 
has multiple self-orthogonal vectors other than the zero 
vector. 

Define the 'conjugate' dual vector for each vector |r) 
in the ortho-nondegenerate basis as 



(r| 



|r> 



|r)-|r) 



(25) 



where it is crucial that |r)-|r) 7^ for (r| to exist. Then, 
the set of dual vectors {(1|, (2|, • • • , (N\} provides a basis 
for the dual vector space V(N,p n )* such that (r\s) = 6 rs . 
Thus, we obtain the set 



{{(l|,(2|,...,(7V|},{|l),|2),...,|iV)}} 
which constitutes a biorthogonal system. 



(26) 



D. Observables 

Given a biorthogonal system, we can define the analog 
of hcrmitian operators via 

JV 

A = J2*k\k)(k\ , a k eGF(p). (27) 

fc=i 

Due to the biorthogonality of the system, \k) is the eigen- 
vector of A with eigenvalue a k . Note that the eigenvalues 
a k are chosen to be elements of GF(p), not GF(p 2 ), i.e. 
they are 'real.' 

In the defining biorthogonal system, the matrix rep- 
resentation of A is diagonal. In a different biorthogonal 
system, say {{<1'|, <2'|, • • • , {N'\}, {|1'>, |2'>, • • • , \N>)}}, 
its matrix representation is 

TV 



(r'\A\s') = J2a k (r'\k)(k\s') 



fc=i 

N 

k=i 



(\r')-\k))(\k)-\ S ')) 
~ k {\r')-\r')){\k)-\k)) ' 



(28) 



which in general is not a hermitian matrix. However, the 
diagonal elements (r'|A|r') are nevertheless 'real' since 
\r')-\k) and \k)-\r') are 'complex conjugates' of each other, 
while \r') ■ \r') and are 'real.' We identify these 

pseudo-hermitian operators with physical observables. 



Since we wish to use Eq. (9) to define the expectation 
value for the observable A, every physical state \4>) must 
have a conjugate dual which we define via Eq. (25). 
Thus, we demand that all physical states belong to some 
biorthogonal system. Essentially, all vectors that are not 
self-orthogonal belong to some biorthogonal system, so 
this requirement is equivalent to dropping all vectors that 
are self-orthogonal from the set of physical states. 

Note that if we multiply \ip) with a scalar, that is, a 
non-zero element of GF(p n ), then its conjugate (tp\ will 
be multiplied by the inverse of that scalar. This will 
leave A and (ip\A\il)) invariant. Thus, we can identify 
all vectors that differ with each other by a multiplicative 
scalar as representing the same physical state, that is, 
all non-zero elements of GF(p n ) can be considered to 
be 'phases.' For V(N,p n ), this means that the set of 
physical states is the non-self-orthogonal subset of the 
projective space 



PG(N-l,p n ) = V(N,p n )\{0} / GF( P n )\{0} 



(29) 



F. Expectation Values 



With the above definitions of observables and physi- 
cal states, we can now calculate the quantity ( - 0|^4|V') € 
GF(p) = Z p for observable A and state \ip). We would 
like to interpret this quantity as the expectation value of 
the observable A. However, if A is to represent a phys- 
ical quantity such as spin, one must map the resulting 
number in GF(p) to a number in R. 

We demand that this map from GF(p) to R be prod- 
uct preserving for reasons that will become clear in the 
following. It is easy to see that the absolute value func- 
tion given in Eq. (14) is a product preserving map for 
any p. For the p = 3 mod 4 case, however, in addition 
to the absolute value function, there is another prod- 
uct preserving map which can be constructed as follows. 
First, denote the generator of the multiplicative group 
GF(p)\{0} by g and express the non-zero elements of 
GF(p) as {p,g 2 7s 3 , • • • = I}- Define: 





+ 1 
-1 



if x = , 
if x = g CVCD , 
if x = g odd . 



(30) 



It is straightforward to show that <p (ab) = ip (a) ip (b) . 

Note that p — 3 mod 4 implies (p — 1) = even and 
(p - l)/2 = odd. Therefore, 



+1 , 



?(-D = <p(g (p - 1)/2 ) = -1, (3i) 

where —1 denotes the additive inverse of 1 in GF(p). 
That is, this function respectively maps — 1, 0, and 1 in 
GF(p) to -1, 0, and 1 in R. 
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We will use this map to give meaning to Eq. (9) as 
an expectation value in the new version of quantum me- 
chanics: 



E{A\i>) = tp 



(<V|i|V> 



(32) 



Using this identification as a starting point in modifying 
ordinary quantum mechanics is a viable alternative to 
specifying a prescription for calculating individual prob- 
abilities for outcomes of measurements, as we mentioned 
earlier. The uniqueness of this map is demonstrated in 
Appendix B. The rule allows us to calculate single and 
joint probability distributions over ensembles. 

An immediate consequence of this rule is noteworthy. 
The uncertainty in the measurement of A will be given 

by 



(AA) 2 = E{A 2 \iP) 
(<^|i 2 



[mm 2 

.))-[^((^|%))] 2 . (33) 



When 
find 



is an eigenvector of A with eigenvalue a, we 



(AA) 2 = V {a 2 )-[ip(a)} 2 = 



(34) 



due to the fact that ip is a product preserving map. Thus, 
if a measurement of an observable is performed on one of 
its eigenstates, the outcome will always be the y-map of 
the eigenvalue associated with that state. If ip were not 
product preserving, this property would not have been 
maintained. 



III. EXAMPLES 



Let us now look at a few concrete examples. 



A. 2D Vector Space over GF(3) 

Consider the 2D vector space V(2, 3) over GF(3) = 
Z 3 = Z/3Z = {0,1,-1}, where we denote the additive 
inverse of 1 as —1 instead of 2. There are 3 2 — 1 = 8 
non-zero vectors in this space which are 



" 1 " 




, \b) = 


' 0" 

1 


, \c) = 


' 1 " 

1 


, !<*> = 


i " 

-I 



\a) 



and their multiples by the 'phase' —1. We find: 



(35) 



\a)-\a) = \b)-\b) = 1, 
\c)-\c) = \d)-\d) = -1. 



Thus, none of the vectors are self-orthogonal, and their 
conjugates are 



(a\ 





1 


Ol 


, (c\ = 


-1 


-1 




Ol 


1 ' 


, (d\ = 


' -1 


1" 



(37) 





fi 


ACTI 


03 


Acr 3 


\a) 





1 


1 





\b) 





1 


-1 





|c> 


1 








1 


|d) 


-1 








1 



TABLE I: Expectation values and uncertainties of spin-like 
observables in biorthogonal quantum mechanics on V(2, 3). 



There are two biorthogonal systems in V(2, 3)* x V(2, 3), 
namely 

{{(a\,(b\},{\a),\b)}} and {{(c|, (d\}, {|c>, \d)}} , (38) 

up to different orderings of the vectors and dual- vectors, 
and signs. All four inequivalent vectors belong to one of 
these biorthogonal systems so they all represent physical 
states. 

We can now construct spin-like observables with eigen- 
values ±1. Since V(2, 3) has only two biorthogonal sys- 
tems, the two possible observables are 



l\a)(a\ -l\b)(b\ 
l\c)(c\ -l\d)(d\ 



1 
-1 

1 

1 



^3 



(39) 



up to signs. By construction, \a) and \b) are respectively 
eigenvectors of a 3 with eigenvalues ±1. Thus, a mea- 
surement of (73 on I a) will always yield +1, while that 
on I b) will always yield —1. Similarly, |c) and \d) are 
respectively eigenvectors of a\ with eigenvalues ±1, so a 
measurement of u\ on |c) will always yield +1, while that 
on \d) will always yield —1. 

On the other hand, the expectation values of &\ and 
o\ for the state \a) are 



Eiaxla) = ^«o|ai|o)) = ^(0) = 0, 
E{a 2 \a) = ip{(a\ai\a)) = V (1) = 1 , 



so 



[A CTl (a)] ; 



E(a 2 \a) 



[E(<yi\a)\ 



(40) 



(41) 



From these expectation values, we can infer the probabil- 
ities of obtaining the outcomes ±1 when o\ is measured 
on \a). Denoting these probabilities as P(±l\a), we must 
have 



P(+l\a) 
P(+l\a) 



l\a) 
l\a) 



(36) which yields 



P(+l\a) = P(-l\a) - 



(42) 



(43) 



Therefore, the measurement of b\ on \a) will yield the 
two outcomes +1 and —1 with equal probability, consis- 
tent with our earlier results in [1] and [2]. Similarly for 
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Acti 


£>~2 


A<7 2 


0"3 


At7 3 


\a) 





1 





1 


1 





\b) 





1 





1 


-1 





c) 


1 








1 





1 


|d) 


-1 








1 





1 


|e> 





1 


1 








1 


l/> 





1 


-1 








1 



TABLE II: Expectation values and uncertainties of spin-like 
observables in biorthogonal quantum mechanics on V(2,9). 



the measurement of &i on \b), and those of 173 on |c) or 
\d). The expectation values and uncertainties for both 
observables and all states are listed in Table I. 

Note that our formalism predicts expectation values 
but do not specify the probabilities directly. The proba- 
bilities must be inferred from the expectations values as 
shown above. Indeed, though we can write 



(a|<7i|a) = (a\(l\c)(c\ 
= l(a\c}(c\a) 



i\d)(d\)\a) 
l(a\d)(d\a) 



(44) 



we cannot associate (a\c)(c\a) = (a\d)(d\a) = — 1 with 
the probabilities of the outcomes ±1. Furthermore, we 
will see in the following that in some cases, the probabili- 
ties cannot be uniquely determined from the expectation 
values. We will argue later that a theory which predicts 
expectation values but leaves the probabilities indetermi- 
nate still makes perfect physical sense. 



B. 2D Vector Space over GF(9) 

Next, consider the 2D vector space V(2, 9) over 
GF(9) = Z 3 [i]. This field consists of 3 2 = 9 elements 
given by {0, 1, — 1, i, —i, 1 + i, 1 — i, — 1 + i, — 1 — i}. 

There are 9 2 — 1 = 80 non-zero vectors in V(2, 9). 
These are the scalar multiples of 80/8 = 10 vectors con- 
sisting of the four listed in Eq. (35) and the following 
six: 



|e) = 
l/) = 



1 

I 

1 

—i 



1.9) = 
\h) = 



1 

i + i 
1 

1 - i 



l*> = 



1 

-I + i 
1 



(45) 

The dot products of these six vectors with themselves are 



|e>-|e) = I/)- 1/) = -I, 

\g).\g) = \h)-\h) = \i)-\i) = \j)-\j) = 



(46) 



As we can see \g), \h), \i), and \j) are all self-orthogonal. 
The conjugates of |e) and |/) are 

<e| = [-1 i] , (/I = [-i-i] ■ (47) 



Thus, in addition to the two biorthogonal systems listed 
in Eq. (38), V(2,9)* x V(2,9) has a third given by 



{{(e|,(/|},{|e),|/)}} 



(48) 



and |e) and |/) are added to the list of physical states. 

The above biorthogonal system contributes a third op- 
erator to the list of spin-like observables in Eq. (39): 



I|e)(e| -ll/X/l = 



= a 2 ■ 



(49) 



By construction, |e) and |/) are respectively eigenvectors 
of er 2 with eigenvalues ±1. The expectation values and 
uncertainties of all three observables for all six states are 
listed in Table II. 



IV. SPIN CORRELATIONS 

In the examples considered above, spin-like observables 
were represented by Pauli matrices, with elements in 
GP(3"), acting on the 2D vector spaces V(2, 3"), n = 1 
or 2. If we associate this model with the spin of one 
particle, two particle spin-states will be represented by 
vectors in V(2,3 n ) ® V(2,3 n ) = V(A,3 n ), n = 1 or 2, 
while the product spins will be represented by Kronecker 
products of the Pauli matrices. In this section, we will 
look at the correlations of these spins. 



A. n = 1 case 

The space V(4, 3) has 3 4 — 1 = 80 non-zero vectors, 
every two of which differ by only a multiplicative phase, 
namely —1, leaving 80/2 = 40 inequivalent vectors. Of 
these, 4 2 = 16 are products of physical states in V(2, 3), 
all of which are also physical in V{A, 3) since 



if {tp\tp) — (<fi\<f>) = i- Of the remaining 
vectors, 16 are self-orthogonal, e.g. 



= I , (50) 
40 - 16 = 24 



" 1" 




" 1 


1 




1 


1 




1 


01 




01 



= 1 + 1 + 1 + Q = , 



(51) 



leaving 24 — 16 — 8 physical entangled states. They are: 

\S) = [014 0] T , 

\(ab)) = [ 1 -l] T , 

\(cd)) = [0 1 1 0] T , 

\(ab){cd)) = [ 1 1 ] T , 

\(ad)(bc)) = [l 1 1 -l] T , 
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(oc)(6d)) = [-1111]' 
|(ac&d)) = [l -1 1 l] 
\(adbc)) = [l 1 -1 l] 



(52) 



where the labeling is based on the transformation prop- 
erty of each state under the group of allowed basis trans- 
formations PO(2, 3). (See Appendix A.l for details.) 

Product spins are represented by <7j <g> &j, i,j = l or 3. 
For product states, the expectation value of product spins 
factorizes due to the product preserving property of <p: 



E{a i a j \i}(j)) = ip 



= £7(^1^) £7(^-10) . 



(53) 



This factorization is necessary if we are to have isolated 
one particle states. Again, the product preserving map 
<p plays a fundamental role. The explicit representations 
of the product spin operators are 



cr 3 (8) cr 3 = 



(54) 



Using these expressions, we can calculate the spin corre- 
lations of this system. 

Let us look at what the Clauser-Hornc-Shimony-Holt 
(CHSH) bound [6] would be. The CHSH bound is the 
upper bound of the absolute value of the following com- 
bination of correlators: 

C(A,a;B,b\t>) 
= E(AB\V) + E(Ab\V) + E(aB\V) - E(ab\V) , (55) 

where A and a are two observables of particle 1, and B 
and b are two observables of particle 2. All four observ- 
ables are assumed to take on only the values ±1 upon 
measurement. For classical hidden variable theory, the 
bound on \C(A, a;B 7 b \^f)\ is 2, while for canonical QM 
it is 2y/2 [10]. 

In the current case, each of the four observables A, 
a, B, and b is either a\ or 0-3. The cases in which the 



operators are the negatives of either a\ or 03 need not 
be considered since 

C(A,a;B,b\V) 
= C(A,-a;b,B\^) = -C(-A,a;b,B\V) 
= c{a,A;B,-b\V>) = -C(a, A ; -B, b . (56) 



To compress our notation, let us define 

Ciju{^) = C{(Ji,(jj\<T k ,o t \if!) 



(57) 



In the current case, there only four possible combina- 
tions of indices: C1313, C1331, C3113, and C3131. Only 
the CHSH correlators for entangled states are of interest, 
since those for the product states cannot exceed the clas- 
sical bound. Furthermore, all eight entangled states can 
be transformed into the singlet state \S) by an appropri- 
ate local PO(2, 3) transformation so one only needs to 
consider correlations for this one state. It is straightfor- 
ward to show that 



From this, wc find 



(S\&3®& 3 \S) 

(5|a 3 ®<Ti|5) 



-I, 
. 



Cl313(S) = C 3 13l(S) = 0, 
Cl33l(S) = C , 3H3(S') = —2 . 



(58) 



(59) 



Thus, the CHSH bound for this model is the classical 2. 

In previous publications [1, 2] we argued that the 
CHSH bound of 2 does not necessarily imply that the 
predictions of the model can be mimicked by a classi- 
cal hidden variable theory. In the current case, however, 
they can be. Let us denote the classical values of o\ and 
03 of particle 1 as X\ and Z\ , and those of the particle 
2 as Xi and Z2, respectively. The first line of Eq. (58) 
implies that the pairs (Xi,X 2 ) and (Zi,Z 2 ) are com- 
pletely anti-correlated. Therefore, the only classical con- 
figurations possible are (X\, Z\\ X 2 , Z 2 ) = (+,+;—,—), 
(+,-;-,+), (-,+;+,-), and (-,-;+,+). To repro- 
duce the second line of Eq. (58), we only need to demand 
that the probabilities of these configurations satisfy: 



= P(+,+;-, -) + P(- -;+,+) 
= P(+,-;-, +) + P(- +;+,-) • 



(60) 



Thus, an entire class of hidden variable mimics exists. 



B. n = 2 case 

The space V(4, 9) has 9 4 — 1 = 6560 non-zero vectors, 
every eight of which differ by only a multiplicative phase, 
i.e. an element of GF(9)\{0}, leaving 6560/8 = 820 in- 
equivalent states. Of the 10 2 = 100 product states, the 
6 2 = 36 products of physical states in V(2, 9) are also 
physical in 1^(4,9). The remaining 64 product states 
are self-orthogonal and unphysical. Of the 820 — 100 = 
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720 entangled states, 216 are self-orthogonal, leaving 
720 — 216 = 504 physical entangled states. These states 
fall into three classes that transform among themselves 
under local PU(2, 9) transformations with 24, 288, and 
192 elements each, as explained in Appendix A. 2. These 
classes can be represented by the following three states 





" 




I 




I 


\S) = 


1 

-1 

_ 


, \T) = 




1 + i 
1 


, \u) = 




1 

i + i 



(61) 



with the duals 



(S\ = [o-iio], 

(T\ = [l l-i l] , 
(U\ = [101 l-i] . 



(62) 



Thus, we only need to calculate the correlators for these 
states to obtain the CHSH bound. Since there are three 
spin observables &i, &2, and a 3 this time, the number of 
possible CHSH correlators is 6 2 =36. 

Let us first look at the correlators involving only o\ 
and a 3 . The correlations for the state \S) are the same 
as those listed in Eq. (58) and (59). Those for the state 
\T) are 



(T|cti®(Ji|T) 
(T|a 3 ®<Ji|T) 
(T|<7 3 ®<j 3 |T) 

from which wc obtain 



(T|<7l®<7 3 |T) 

I, 
, 



-1 



(63) 



Cl 3 l 3 (T) 

C 3 n 3 {T) — C 3 i3i(T) = 1 , 
Ci33i(T) = -3. 

Similarly, for the state \U) we have 

{U\o x ®oi\U) = {U\a x ®a 3 \U) 

= (U\a 3 ®a 3 \U) = -1, 
<I7|* 3 ®*i|E0 = I, 



and 



C 13 i 3 (U) = C 3U3 (U) 
C 1331 (U) = -4. 



C 3 13l(^) = 



(65) 



(66) 



As can be seen, the absolute value of the correlator C1331 
for the states \T) and \U) exceed not only the classical 
bound of 2 but also the Cirel'son bound of 2\pl. In a 
similar fashion, we have scanned all 36 spin combinations 
for the three states and have obtained the tally shown in 
Table III. Thus, wc find that the CHSH bound for this 
model is 4. 

Unlike the n = 1 case, which had a CHSH bound of 2, 
the above correlations cannot be reproduced by any clas- 
sical hidden variable theory. For instance, the first line of 
Eq. (63) demands that the pairs (Xi,X 2 ) and [X\,Z-i) 



state 





1 


2 


3 


4 


IS) 


6 


24 


6 








|T> 


6 


18 


6 


6 





|C0 


12 


12 


4 


4 


4 



TABLE III: The number of CHSH correlators with the re- 
spective absolute values for the three states \S), |T), and \U). 



are completely anti-correlated, while the second line de- 
mands that the pair (Zi,X 2 ) is completely correlated. 
But then X\ = ±1 would imply X2 = +1 and Z2 = TL 
the first of which implies Z\ = =pl. Therefore, the pair 
(Zi,Z 2 ) must also be completely correlated which con- 
tradicts the third line of Eq. (63). Similarly, Eq. (65) de- 
mands that the pairs (X\, X 2 ), (Xi, Z2), and (Zi, Z 2 ) are 
completely anti-correlated, while (Zi,A" 2 ) is completely 
correlated. But then Xi = ±1 would imply X 2 = Tl and 
Z 2 = TL the former of which implies Z\=^\ while the 
latter Z\ = ±1, leading to a contradiction. Of course, 
this is not surprising since the CHSH bound for classi- 
cal hidden variable theories is 2. The unexpected result 
is that the CHSH bound of our model also exceeds the 
quantum Cirel'son bound of 2^/2. In the next section, 
we will take a careful look at how this comes about. 



V. EXPECTATION VALUES WITHOUT 
DEFINITE PROBABILITIES 

In canonical QM, the states that correspond to \S), 
\T), and \U) are 



(64) \S) 



1 

T2 








1 




1 


1 

-1 


' \^=\ 




l + l 





1 







1 







(67) 

Calculating the correlations of canonical spin <7j for the 
state \S) in canonical QM, we find 

(S|oti®cti|S) = (S\& 3 ®a 3 \S) = -1 , 
(S|ori<8>* 3 |5) = (S|? 3 ®*i|5) = 0, (68) 

which agree with those for 15*) in Eq. (58) via the product 
preserving map ip. For \T) and \U), however, we find: 



(T|ai®CTi|f) 

(T|ct 3 ®cti|T) 
(f\a 3 ®a 3 \f) 



(T|CTi®a 3 |T) = 
1 

_ 2 ' 



= 0, 



(U\a 3 ®a 3 \U) 
1 

_ 2 ' 



(69) 
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Thus, the correspondence here is 

-I \ , I ^ ~\, (70) 

which is to be expected since 1 2 = 2 = — 1 in GF(3). 
So the large correlation is due to the fact that GF(3) 
has only three elements {—1,0,1} which are mapped to 
{—1,0,1} € M. by the product preserving map </?. The 
fact that the only spin-correlations possible are or ±1 
will of course persist for larger values of p — 3 mod 4 as 
long as we use ip. 

What are the corresponding probabilities? Let us take 
the spins in the Z-direction, 173 ® 03, as an example. 
The probabilities of the outcomes (03(73) = (++), (H — ), 

( — h), and ( ) in canonical QM are listed in Table IV. 

As can be seen, they reproduce the correlations listed 
above as they should. 

In our 'mutant' biorthogonal quantum mechanics, 
however, the probabilities of individual outcomes are ill 
defined as discussed above. Taking the point of view that 
the probabilities must be inferred from the expectation 
values, we have the constraints 

P{+ + |T) + P(+ - \T) + P{- + \T) + P(- - |T) = 1 , 
P(+ + \T) - P(+ - \T) - P{- + \T) + P{- - \T) = , 

(71) 

for |T), and 

P(+ + \U) + P(+ - \U) + P(- + \U) + P{- - \U) = 1 , 
P(+ + \U) - P(+ - \U) - P(- + \U) + P(- - \U) = -1 , 

(72) 

for I U) . These constraints imply 

\ = P(+ + \T) + P(--\T) 
= P(+ - \T) + P(- + \T) , 

= P(+ + \U) + P(--\U) , 

1 = P(+-\U) + P(- + \U) , (73) 

but beyond this the probabilities cannot be specified. 
Therefore, though our formalism predicts definite expec- 
tation values, it leaves probabilities indeterminate. Phys- 
ically, we interpret this to mean that if the same mea- 
surement is repeated many times, the average of the out- 
comes will converge to the predicted expectation value, 





++ 


+- 


-+ 




E.V. 


\S) 





1 

2 


1 

2 





-1 


\f) 


1 

4 





1 

2 


1 

4 





P) 


1 

4 





1 

4 


1 

2 


1 

+ 2 



TABLE IV: The probabilities of the four possible outcomes 

++, H — , — h, and in canonical quantum mechanic when 

03 (g><5"3 is measured on the canonical states \S), \T), and \ U). 



while the frequencies of each outcome will continue to 
fluctuate. 

This indeterminacy is characteristic of the approach 
used here, and can be understood more generally by re- 
examining the defining relation between expectation val- 
ues and probability distributions. In conventional QM, 
it is possible to construct the probability distribution for 
the measurement outcomes of some observable through 
use of the system of equations formed by the expectation 
values of the powers of the observable in question. This 
is not possible for spin observables in the model under 
consideration due to the cyclic nature of the underlying 
field. Explicitly, the system of equations: 

E{A\$) - 5>P(c#), 
E(A 2 \iP) = f> 2 P(c#), 

a 

E{A N \%jj) = J2a N P(a\<P), (74) 

a 

will be singular if V is greater than the least common 
multiple of the multiplicative orders of the eigenvalues 
a of A since the cyclic nature of the field is necessarily 
shared by the eigenvalues when the product preserving 
map also preserves the eigenvalues. In our examples, us- 
ing GF(3) as the 'real' field, the eigenvalues of spin ob- 
servables, {+1, —1}, have multiplicative orders no greater 
than 2. Thus, when we form a four level system by en- 
tangling two particles, we find that the system of equa- 
tions needed to solve for the probabilities of these four 
measurement outcomes is singular and cannot be used to 
assign consistent probabilities. 

In Ref. [3], we conjectured that a 'doubly' quantized 
theory may predict super-quantum correlations with a 
CHSH bound which exceeds the Cirel'son value of 2\/2. 
A state in such a theory can be thought of as a 'su- 
perposition' of various 'singly' quantized states, each of 
which predicts definite probabilities. A 'measurement' 
in a 'doubly' quantized theory can be expected to col- 
lapse the 'doubly' quantized state to a 'singly' quantized 
one, selecting a particular probability distribution from 
all possible ones. Every 'measurement' will lead to a 
different probability distribution, so no definite probabil- 
ity will be predicted. These considerations suggest that 
biorthogonal QM is a candidate model for such a 'doubly' 
quantized theory. 

VI. DISCUSSION 

One of the simplest realization of how quantum the- 
ory differs from its classical counterpart is given by the 
celebrated Bell inequalities, or its slightly generalized ver- 
sion, the CHSH inequalities [6-10]. According to these 
inequalities the classical and quantum physics are clearly 
separated by O(l) effects. It has been pointed out in the 
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literature that the purely statistical reasoning leads to 
the maximal "super-quantum bound" of 4 [13]. In one of 
our previous papers we have pointed out the special na- 
ture of such a super-quantum theory [3]. Given the fact 
that the CHSH inequalities rely on the knowledge of ex- 
pectation values (and not probabilities) in this paper we 
have focused on the requirement that expectation values 
of a super-quantum theory should satisfy the bound of 4. 

Note that our present work is distinguished from other 
efforts that try to eliminate theories which violate the 
quantum bound or which claim the uniqueness of the 
canonical complex quantum theory because of the sup- 
posed unphysical nature of super-quantum theories (see 
[14]). As is well known, the expectation values and the 
probabilities are related by a quadratic map in canonical 
quantum theories and its real counterparts [15]. That 
this map is quadratic can be argued on general grounds, 
and the robustness of the Born rule [16], by pointing out 
the generic nature of the Fisher metric on the space of 
measured events [17]. 

The CHSH observable relies only on the computation 
of the expectation values. In order to achieve the super- 
quantum bound of 4, one immediately realizes (at least 
on a heuristic level) that the expectation values should 
be "mutated" so that the last term in the CHSH observ- 
able changes its sign. Given the canonical relation be- 
tween the expectation values and the probabilities, such 
a "mutation" of the computation of the expectation val- 
ues would, at least naively, influence the probabilities as 
well. This is precisely what we find in a concrete mathe- 
matical model explored in this paper: the CHSH observ- 
able computed in the mutant quantum mechanics over 
the finite field GF(9) is explicitly equal to 4, which in 
turn implies that the probabilities are indeterminate in 
such a super-quantum theory. 

Indeterminate probabilities are a consequence of our 
construction and, in this particular case, a necessary fea- 
ture of such a super-quantum theory. Note that this 
statement also goes against some efforts in the founda- 
tions of quantum theory, which try to base the canonical 
complex quantum theory solely on the concept of prob- 
ability (see for example [18]). Our point is that even 
though canonical quantum theory might be solely based 
on the concept of probability, super-quantum theory does 
not have to be. This reinforces the experience of mod- 
ern QFT (especially the conformal QFT's) in which one 
operates only with correlation functions. 

In Appendix A, we show that, in the context of Galois 
biorthogonal QM, the projective orthogonal and the pro- 
jective unitary groups play the natural role of the orthog- 
onal and unitary groups of canonical QM. This maintains 
a parallel with our previous papers on Galois field QM 
[1, 2] where we have shown that the complex (and real) 
projective spaces, which define the geometry of canonical 
quantum theory, can be naturally replaced by their finite 
projective counterparts. Similarly, in this work, the or- 
thogonal and unitary groups that define the invariancc 
of expectation values in the real and complex quantum 



theories are replaced by their projective counterparts. It 
is of course tempting to contemplate that the general 
structure of biorthogonal systems, the graded valuation 
of expectation values, and the indeterminate nature of 
probabilities is valid for more general constructions of 
super-quantum theories, including the ones that we ex- 
pect to be relevant in quantum theory of gravity. 

To summarize: in this paper we have presented per- 
haps the simplest model for quantum super-correlations. 
Quantum super-correlations are realized in the model to- 
gether with a signature feature: the physics of the model 
is entirely determined in terms of expectation values, 
whereas the probabilities are, in general, indeterminate. 
This feature is actually quite natural (and desirable) from 
various point of view suggested by different modern av- 
enues of fundamental physics. 

We note that the fact that the probabilities are inde- 
terminate in our explicit construction also meshes well 
with some expectations from various attempts at quan- 
tum theory of gravity (including the ones in which con- 
formal field theories are used to define a quantum theory 
of gravity in particular asymptotic geometries.) Indeed, 
that fundamental quantum theories can be defined in 
terms of expectation values (which is most obvious in the 
path integral formulation), is a feature found in modern 
conformal field theories, which are quantum field theory 
formulated from a purely algebraic viewpoint, without 
the use of Lagrangians (or Hamiltonians) or Feynman 
rules. For example, the familiar S-matrix of the canoni- 
cal quantum field theory, which comes about from com- 
pounding expectation values (correlation functions) with 
wave-functions of external probes, is not a well-defined 
concept in conformal field theory. As is well known, con- 
formal field theories, can be dual to (quantum) gravi- 
tational theories in certain background (the AdS spaces 
[19], and also in the context of the observed cosmological 
de Sitter spacetimes [20]). 

Thus, this feature should be relevant in the context 
of quantum gravity as well. Indeed, different approaches 
to non-perturbative quantum gravity and quantum cos- 
mology [21-23], suggest that the individual probability 
for specific measurements could be indeterminate, and 
that the observables in that context are different from 
the usual observables found in the canonical quantum 
theory. The model considered here should be viewed as 
a concrete realization of this general expectation. 

The model sheds new light on the foundations of quan- 
tum theory, and attempts to understand the simplest set 
of reasonable axioms that lead to canonical quantum the- 
ory, which could lead to natural generalizations of quan- 
tum theory expected in the context of quantum theory 
of gravity [23, 24]. 

Finally, we note that this work presents an alternative 
pathway to constructing a quantum theory on a vector 
space without an inner product from the one introduced 
in Refs. [1,2]. Application of the two constructions to Ba- 
nach spaces [25] would be a natural place to further clar- 
ify the difference between the two approaches, do away 
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with the product preserving map from GF(p) to R, and 
search for models which may serve as closer representa- 
tions of reality where various quantum gravitational ideas 
discussed above can be explored. 

We will return to these, and related issues in future 
works. 
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Appendix A: Group of Basis Transformations 
1. V(2,3) case 

There are only two biorthogonal systems in V(2, 3)* x 
1/(2,3) listed in Eq. (38), up to ordering of the vectors 
and multiplicative phases. Thus, the allowed bases of 
1/(2,3) are 



±{\a),±\b)}, 
±{\c),±\d)}, 



±{|6),±|a)}. 
±{\d),±\c)}. 



(Al) 



Thus, the group of all possible basis transformations con- 
sist of sixteen matrices given by 



e «-> ± 

(cd) <-> ± 

{ac){bd) ± 

(acbd) <H- ± 



10 , 1 

- , lab) -H> ± ~ 

1 10 

- , (ab)(cd) «-» ± ~ 

0- 1 10 

1 1 , (ad)(bc) o ± _i 1 

1- I II 

- , (adbc) ± _ 
1 1 J V ' [-II 

' (A2) 
However, since we identify vectors that only differ by mul- 
tiplicative phases as representing the same physical state, 
we identify the matrices that only differ by a multiplica- 
tive phase as representing the same transformation on 
the projective space PG(1, 3), each of which corresponds 
to a permutation of the vector labels a, b, c, and d as 
indicated above. These eight transformations constitute 
the projective orthogonal group PO(2,3) = D 4 , namely, 
the group of 2 x 2 matrices O with elements in GF(3) 
which satisfy the condition 




FIG. 1: The correspondence between the dihedral group D4 
and the projective orthogonal group PO(2,3). Every D4 ro- 
tation of the quadrangle corresponds to a permutation of the 
four vertex labels abed belonging to PO(2, 3). 



with matrices which differ by a sign identified. This 
group is a subgroup of the projective general linear group 
PGL(2, 3) ^ £4. 

The isomorphism between PO(2,3) and D4 is imple- 
mented by labeling the four corners of a square as shown 
in Fig. 1. Every rotation of the quadrangle in D4 leads 
to a permutation of the four vertex labels, which is the 
corresponding element of PO(2, 3). The two spin observ- 
ables (Ti and (T3 transform under PO(2,3) permutations 



(06) 
(cd) 

(ab) (cd) 

(ac) (bd) 

(ad) (bc) 
(acbd) 
(adbc) 



01 -> 01 j 03 ^ 03 3 
01 — > 01 , 03 — > —03 i 

01 — > — 01 ) 03 -> 03 5 
01 —01 , 03 — > ^03 

01 ->• 03 > 03 01 i 
01 — > —03 ) 03 — > —01 

01 -> -03 ) 03 -> 01 ! 
01 — >• 03 ! 03 — >• -01 ! 



(A4) 



just as they should under rotations of the quadrangle. 

The eight elements of PO(2, 3) fall into five conjugacy 
classes given by 

{e}, 

{(ab)(cd)} , 
{(06), (cd)}, 

{(ac)(6d), (ad) (6c)} , and 

{(ac6d), (ad&c)} . (A5) 

The eight physical entangled states in V(2, 3) x V(2, 3) = 
V(4, 3) also fall into five classes that transform among 
themselves under global PO(2, 3). They can be classified 
and labeled according to their transformation properties 
under the full global PGL(2, 3). 



\S) 



[0 1 -1 o]' 



o T o 



z " = l2x2 



(A3) 



\(ab)) = [100 -I] 
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\(cd)) = 


[o 


1 


1 





] T . 


ab){cd)) = 


[ 1 








1 


r. 


ad)(bc)) = 


[ 1 


1 


1 




i] t 


ac){bd)) = 




1 


1 


1 


i] t 


\(acbd)) = 


[ 1 




1 


1 


i] t 


\(adbc)) = 


[ 1 


1 




1 





mation are represented by the following matrices: 



(A6) 



Here, \S) is the singlet state which is invariant under 
all transformations in PGL(2,3). The state \(ab)(cd)) 
is also a singlet under PO(2, 3) transformations, but 
transforms into \(ac)(bd)) and \(ad)(cd)} under the full 
PGL(2,3). The other states transform in pairs under 
PO(2,3), falling into the same classes as the PO(2, 3) 
transformations themselves as listed in Eq. (A5). 

Under local PO(2, 3) transformations, that is, 
PO(2, 3) transformations acting on only one of the 
V(2,3) vector spaces in U(4,3) = V(2,3) x V(2,3), all 
eight states fall into the same class and can be trans- 
formed into the singlet state \S). Explicitly, we have: 

\S) = (cd)i|(cd)) = (cd) 2 \(cd)) 

= (ab^Kab)) = (ab) 2 \(ab)) 

= (a6)i(cd)i|(a6)(cd)) - (ab) 2 (cd) 2 \(ab)(cd)} 

= (ac)i(6d)i|(oc)(6d)> = (ac) 2 (bd) 2 \(ac)(bd)) 

= (ad)i(6c)i|(ad)(6c)> = (ad) 2 (bc) 2 \(ad)(bd)) 

= (acbd)i\(acbd)) — (acbd) 2 \(adbc)} 

= (adbc)i\(adbc)) = (adbc) 2 \(acbd)) , (A7) 

where the subscript indicates which V(2, 3) space the 
transformations are acting on. 

The above considerations indicate that it suffices to 
calculate the CHSH bound for only the singlet state \S). 



2. V(2,9) case 

In V(2,9), we have three biorthogonal systems listed 
in Eqs. (38) and (48). Unlike V(2,3), this space has 
unphysical self-orthogonal vectors so some care is nec- 
essary in listing possible bases since the dot product is 
not invariant under generic basis transformations, and 
a non-self-orthogonal vector may be mapped to a self- 
orthogonal one. Using the notation of Eqs. (35) and (45), 
the allowed bases are 



r?{|a),±|6)}, 
r?{|6),±|a)}, 
V{\c),±\d)}, 
V{\d),±\c)}, 
V{\e),±\f)}, 
V{\f),±\e)}, 



V{\a),±i\b)}, 
V{\b),±i\a)}, 
V{\c),±i\d)}, 
V{\d),±i\c)}, 
r?{|e),±i|/)}, 
V{\f),±i\e)}, 



(A8) 



e -H- r] 
(cd)(ef) f+r? 
(acbd) o i] 
(adbc) «-> 77 
(aebf) «-» i] 
(afbe) f] 
(cedf) o i] 
(cfde) <-> 77 

(ace) (bdf) «-» 77 

(acf) (bde) «-» 77 

(aec) (bfd) <h> r\ 

(aed) (bfc) «-» 77 



- - , (ab)(ef) O 77 - 

1 1Q 

" , (a&)(cd) <-> 77 - - 

~ ~- . {ac)(bd)(ef) ^77 - _- 

_^ J , (ad)(bc)(ef) <-> 77 ~j J 

1 * , (ae)(bf)(cd) O 77 1 ~- 

2 1 7 —1 

1 ~" , {af)(be){cd) <-> 77 1 1 

-7 1 -7-1 

q ^ . (a&)(ce)(d/) o 77 ° ^ 

n 2 , (a6)(c/)(de) o 77 2 * 
0—7 10 

7 ~- , (adf)(bce) Or/ ; 1 
; - , (ade)(bcf) <-> 77 ~; 1 

1-7 1 7 

- - , (afd)(bec) «■ 77 ~- - 
7—7 11 

- ~ , (afc)(bed) o 77 - - 
77 — 1 i 

(A9) 

Identifying matrices that differ by a multiplicative phase, 
we obtain a group of basis transformation with 24 el- 
ements, each of which corresponds to a permutation of 
the vector labels abedef as indicated above. This group is 
the projective unitary group PU(2, 9) consisting of 2 x 2 
matrices J7 with elements in GF(9) which satisfy the con- 
dition 



U^U = ±1, 



(A10) 



where 77 is an arbitrary phase, that is, an element of 
GF(9)\{0}. Thus, the group of allowed basis transfor- 



with matrices that differ by a multiplicative phase iden- 
tified. This group is a subgroup of PGL(2, 9) which is 
isomorphic to the octahedral group O, which is also iso- 
morphic to S4. 

The isomorphism between PU(2, 9) and the octahedral 
group O is implemented by labeling the six vertices of the 
octahedron as shown in Fig. 2. Every rotation of the oc- 
tahedron in O will lead to a permutation of the vertex 
labels corresponding to an element of PC/(2,9). For in- 
stance, the 180° rotation around the x-axis lead to the 
permutation (ab)(ef). The spin observables transform 
under PU(2, 9) as 

e : Oi <7i , <7 2 -> <7 2 , <7 3 -> <7 3 , 

(ab)(ef) : a\ -> a\ , a 2 -> -<r 2 , <r 3 -> -<r 3 , 
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Z 




x 



FIG. 2: The correspondence between the octahedral group O 
and the projective unitary group PU(2,9). Every O rotation 
of the octahedron corresponds to a permutation of the six 
vertex labels abcdef belonging to PU(2,9). 



(ab)(cd) 


0*1 


-O-l 


! 02 - 


"> 02 , 0*3 -> -03 j 


(cd)(ef) 


01 


-5- -<Ti 


> 02 - 


-> -^2 , 03 -> 03 ; 


(aebf) 


0*1 


-> 0*1 , 


0*2 -> 


— 03 i 03 ^ 0*2 ; 


(afbe) 


0*1 


-> 01 , 


0*2 — > 0*3 , <7 3 — > — <7 2 , 


(acbd) 


0*1 


-> -CT3 


! 02 - 


"> 02 , 0*3 ^ 0*1 ? 


(adbc) 


0*1 


-> 0*3 , 


<7 2 -> 


0*2 , 0*3 -> -0*1 ; 


(cedf) 


0*1 


-> 0*2 , 


(72 -> 


-0*1 , 0*3 -> 0*3 , 


(cfde) 


0*1 


-> — CT2 


! 02 - 


"> 0*1 j 0*3 -> 0*3 , 


(ae)(bf)(cd) 


0*1 


— > — CTl 


i 02 - 


0*3 j 0*3 -> 0*2 , 


(af)(be)(cd) 


0*1 


-> -CTl 


i 02 - 


-> —(73 , (73 — > —(72 


(ac)(bd)(ef) 


0*1 


-> 0*3 , 


0*2 -> 


-0*2 , 0*3 -> 0*1 , 


(ad)(bc)(ef) 


0*1 


-> -CT3 


, (72 -> -(72 , 03 -> 0" 1 


(ab)(ce)(df) 


0*1 


0*2 , 


& 2 -> 


0*1 , 0*3 -> -0*3 , 


(ab)(cf)(de) 


0*1 


-> -02 


! 02 - 


-» —(7i , (73 — > — (73 


(ace)(bdf) 


0*1 


-> 0*2 , 


02 ^ 


0*3 , 0*3 -> 0*1 , 


(adf)(bce) 


0*1 


0*2 , 


0*2 -> 


— 0*3 j 0*3 -> -0*1 j 


(acf)(bde) 


0*1 


-> -CT2 


! 02 - 


"> -0*3 ! 0*3 ->• 0*1 j 


(ade)(bcf) 


0*1 


-> -02 


! 02 - 


0*3 j 0*3 -> —0*1 j 


(aec)(bfd) 


0*1 


-> 0*3 , 


02 -> 


0*1 , 0*3 — ► 0*2 , 


(afc)(bed) 


0*1 


— > (73 , 


CT2 -> 


— 0*1 , 0*3 -> —0*2 , 


(aed)(bfc) 


01 


-> -0"3 


i 02 - 


-> -<7i , <7 3 -> (7 2 , 


(afd)(bed) 


CTl 


-> -CT3 


! 02 - 


"> 0*1 j 0*3 -> -0*2 , 



(All) 



just as they should under the corresponding rotations of 
the octahedron in 3D space. 



The five conjugacy classes of PU(2, 9) = O are 
{e}, 

{(ab)(ef),(ab)(cd),(cd)(ef)} , 

{(acbd), (adbc), (aebf), (afbe), (cedf), (cfde)} , 

{(ac)(bd)(ef), (ad)(bc)(ef), (ae)(bf)(cd), 

(af)(be)(cd),(ab)(ce)(df),(ab)(cf)(de)} , and 
{(ace)(bdf), (adf)(bce), (acf)(bde), (ade)(bcf), 

(aec)(bed), (add) (bed), (add)(bfc), (afc)(bed)} . 

(A12) 

The 504 physical entangled states in V(2,9) x V(2,9) = 
V(A, 9) also fall into classes that transform among them- 
selves under global PU(2, 9) transformations. Since we 
cannot list all 504 states here, we will only mention that 
they fall into 17 classes of 24 elements each, 4 classes of 
12 elements each, 4 classes of 8 elements each, 2 classes 
of 6 elements each, 1 class of 3 elements, and the singlet 
state \S) = [0,1,-1,0] T . 

This can be verified by a direct search, or through 
use of Burnside's lemma and the orbit-stabilizer theo- 
rem. For a group, G, acting on a set, X, a subset that 
is preserved by the action of the entire group is called an 
orbit. The set of orbits forms a partition of the set X. 
To calculate the number of orbits, here denoted \X/G\, 
one can use Burnside's lemma: 

|G| 9 ec 

where X 9 is the set of elements in X that are invariant 
under the action of g. 

For global PU(2, 9) transformations, Burnside's 
lemma indicates that there should be 29 orbits in the 
set of entangled states. To calculate the length of these 
orbits, one could use the orbit-stabilizer theorem, which 
states that the order of the orbit containing an element 
is equal to the order of the group divided by the order 
of the stabilizer subgroup of that element. The stabilizer 
subgroup of an element is the subgroup under which that 
element is invariant. 

This computation indicates that there are 408 states 
that belong to orbits of order 24, 48 states that belong to 
orbits of order 12, 32 states that belong to orbits of order 
8, 12 states that belong to orbits of order 6, 3 states that 
belong to orbits of order 3, and 1 states that belongs to 
an orbit of order 1. 

Under local PU(2, 9) transformations, the same 504 
entangled states fall into three classes with 24, 288, and 
192 elements each. Again, this result can be arrived at 
through a manual search or through the group theoretic 
means mentioned above. Representative elements from 
the three classes can be taken to be: 

\S) = [Q 1 -1 o] T , 

|T) = [10 1 + i l] T , 

1 17) = [l 1 l+i] T • (A14) 

Therefore, to obtain the CHSH bound, we only need to 
calculate the correlators for these three states. 
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Appendix B: Uniqueness of the Product Preserving 

Map 

Is the (p function defined in Eq. (30) the only function 
that allows us to calculate real expectation values? Yes, 
as we will see via the following argument. 

What are the physical requirements on ipl Primarily, 
it must be a map from GF(p) to K, as we assume that 
the results of measurements are real numbers. For the ex- 
pectation value of the identity operator to be 1, we must 
have that ip(l) — 1. Likewise, for the expectation value 
of the zero operator to be 0, we must have that <p(0) = 0. 
When considering two particle states, if we require that 
the expectation values of product states should factor- 
ize, then ip must respect multiplication. This is done if 
the image of GF(p)\{0} is homomorphic to GF(p)\{0}. 
Since GF(p)\{0} is cyclic, any group homomorphic to 
it must also by cyclic. The only cyclic, multiplicative 
subgroups of K are {+1}, {0}, and {+1, —1}. 

As the image must contain +1, {0} is excluded. If 
we choose {+1,-1} as the image of GF(p)\{0} under <p, 



we are allowed to interpret eigenvalues of observables as 
the expectation values of the corresponding eigenstates. 
Thus, the image of GF(p)\{0} should be {+1, -1} and <p 
should be surjective between GF(p)\{0} and {+1,-1}, 
to ensure the presence of —1 as an expectation value. 

For ip to be a such a surjection, it must be true that 
ip(g) = — 1 whenever g is a multiplicative generator of 
GF(p)\{0}. It then follows that all even powers of g 
should map to +1. Thus, the kernel of <p must contain all 
(p — l)/2 even powers of g. Note that it does not matter 
which generator is chosen since any given generator is an 
odd power of each of the other generators. 

The kernel of a group homomorphism, the set of ele- 
ments that map to the identity, is a subgroup. In order 
for ip to be surjective from GF(p)\{0} to {+1, —1}, its 
kernel must be a proper subgroup of GF(p)\{0}. As 
we have shown that the kernel must contain half of the 
elements of GF(p)\{0}, it can only contain those ele- 
ments, as the order of the kernel must divide the order 
of GF(p)\{0}. Therefore, <p as defined in Eq. (30) is the 
only map that fits the relevant criteria. 
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